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ABSTRACT 


A field theory approach is used to develop the volt- 
ampere and force equations of an infinitely long tubular motor 
for quasi-static conditions. A circuit model is developed 
which includes the spatial harmonics of the electromagnetic 


fields. 


Analysis of the effect of stator winding distribution, 
air gap width, rotor sleeve thickness and constant clearance 
shows that the air gap width should be as small as possible and 
the rotor sleeve thickness should be chosen to give optimum 


resistances for the best machine performance. 
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CHAPTER i 


INTRODUCTION 


1.1. Linear Induction Machine 

Extensive pioneering work on linear induction motors was done by 
Laithwaite (2) about a decade ago. Today many engineers are involved 
in projects on transit systems, using single-sided and double-sided 
linear machines as a means of propulsion. 

A linear induction motor can be simply described as a conven- 
tional rotary machine (Fig. 1-1) which has been cut along an axial plane 
and opened out flat (Fig. 1-2) to give a single-sided linear induction 
motor (LIM). Addition of a second stator (Fig. 1-3) leads to a double- 
sided LIM. Considering relative rotor and stator lengths, LIMs may be 
divided into two classes. These classes have been termed as short- 
stator (rotor longer than stator) and short-rotor (stator longer than 
rotor) machines. 

If the model in Fig. 1-2 is re-rolled about an axis parallel to 
the direction of the field motion, the resulting model (Fig. 1-4a) repre- 
sents the basic structure of a tubular motor. The winding of the tubular 


motor consists of an array of coils (Fig. 1-4b), having no end turns. 


STATOR 


FIG. 1-1 CONVENTIONAL INDUCTION MACHINE 
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. 1-4 TUBULAR MOTOR 


(a) basic structure 
(b) array of coils 


As can be seen from Fig. 1-4a, where N and S represent alternate 
N and S poles, the tubular motor is an axial-flux machine. The flux 
is forced axially along the rotor, therefore the rotor must contain 
iron if the magnetic circuit is to be at all reasonable. 

A number of liquid metal induction pumps employing the use of 
tubular stator windings have been devised for pumping liquid bismuth, 
sodium or potassium. Utilizing the same winding arrangements and 
replacing the liquid metal with a copper or aluminium sleeve may result 
in two different methods of constructing the rotor. In general, the 
rotor sleeve may be permanently affixed to the rotor iron or the rotor 
iron may be held by some external means while the rotor sleeve is 
allowed to glide (6). The analysis in the following chapters is 


applicable for both types of rotor structure. 


1.2. Method of Analysis 

The classical analysis of induction machines (11) applied to the 
tubular motor would not be adequate because of considerations such as 
large air gaps and lack of definitive rotor conductors. Therefore, a 
field theory approach (such as used by (1,3)) is necessary and is used 
to develop the force expressions for steady state sinusoidal excitation. 

Energy conversion characteristics of a polyphase induction machine 
with symmetrical impedances can be derived from an equivalent two-phase 
induction machine by a change of variables (5). For this reason a 
two-phase tubular motor is considered for simplicity. 

A lumped parameter circuit model based on the field theory is 


developed as the circuit model offers advantages in the computation of 
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input and output quantities under steady state conditions. 

Tubular motors operate at low frequencies (usually 50 or 60 Hz), 
therefore a quasi-static approach (1,3,5) is used to determine the 
fields in the air gap. Equivalent inductances of the machine are found 
from the magnetic storage energy and then the equations of motion are 
derived. Using linear axis transformations, the spatial dependance of 
the rotor inductances is removed and an equivalent circuit is found. 
Spatial harmonics of the current densities are included, resulting in 
a model which consists of a set of harmonic machines with series- 
connected stator windings and a common output shaft. 

Calculations of the forces are done on a digital computer and 
are compared for different values of air gap, stator winding distri- 


bution and rotor conductivity. 


1.3. Objective of Thesis 


The steady-state characteristics of a tubular motor are derived 
from equations of motion based on the field theory. An APL program is 
set up to calculate steady state force, efficiency, input current and 
harmonic forces for different values of rotor slip. The effect of the 
rotor sleeve thickness, air gap thickness and stator winding distribu- 
tion on the machine performance is examined. This is done for two 


kinds of rotor sleeve material - copper and aluminum. 
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CHAPTER II 


ELECTROMAGNETIC FIELD 


Electromagnetic field equations are developed from Maxwell's equations 


for a described model. In essence, the field equations are partial 


differential equations. These are solved subject to the given boundary 


conditions. 
ape a Assumptions 


The following simplifying assumptions will be made in deriving 


the field equations: 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


The stator and rotor iron cores are laminated so that conductivity 
(co) is zero. They have a permeability u for which u >> uo holds 
and therefore »p is assumed to be infinitely large. 

The tubular motor (Fig. 2-1) is replaced by the model shown in 
Fig. 2-2, and the coordinate system is as shown. This model is 
adequate if the air gap is small with respect to the rotor 
diameter. 

The longitudinal effect of the motor is neglected (both the rotor 
and stator are infinitely long in the z direction). 

The air gap is uniform. 

For the operating frequency of 60 Hz the current density of both 
the rotor and stator is uniform in the z direction, that is, the 
skin effect is neglected. 


The stator conductors are not skewed. 
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FIG. 2-1 TUBULAR MACHINE 
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FIG, 2-2 CARTESIAN CO-ORDINATE SYSTEM 
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2.2. Stator and Rotor Current Densities 

Three different types of stator windings are considered. The 
type I winding (14) as shown in Fig. 2-3 is the simplest and least 
expensive winding to make. The type II and III windings produce the 
current density distribution as shown in Fig. 2-5 and Fig. 2-6 


respectively. 


CURRENT iN | 
DENSITY | 


FIG. 2-3 STATOR WINDING (one phase) TYPE I 


The stator winding is not replaced by a current sheet, thus the 
winding has a finite thickness which is represented by the value of 
ly. = y,| in Fig. 2-2. In general, the stator current density will be 
of the form as shown in Fig. 2-4. 


i 


CURRENT 
DENSITY 


FIG. 2-4 GENERAL STATOR CURRENT DENSITY (one phase) 


I BavT fseadq smo) DHIGVIW SOTAT2 E-S .O1F 


ody auld ,Iseda scisv1v2 8 yd beos{qet Som et guibatw s0IsIe eT 
to eeley S115. yd baanseerqot 2f dotdw aeomlotds stint? 5 ssii gabbalw 


od Iftw wWisnsb dnerzua roinse sf3 ,{nxraang nt .S-S ght nt Pi - exl - 


-d-S . 32% ot awole as myod sis Jo . . 
i 
1 | arenes 


Following the method introduced by Hague and developed by Saunders 


(3) and Nasar (4), the stator current density is defined as 


merece elk aa 2 
J35 a Gy, 5 y,) (amperes/meter~) (2-1) 


th 


The k~ stator current density pulse (Fig. 2-4) is generally 


expressed in the Fourier series 


Nem 
k 
hy (x) = a on 2 Vay COn yx ta Dy sin A,X) 
where 
2TV 
A) ri cm Wy = ¥3 = Yo 
and 
5 Tyr/2 
agpe bie i £ (x) cos A. x dx 
ig 
2. 
CURRENT h 
DENSITY : 


FIG. 2-5 STATOR WINDING (one phase) TYPE II 
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Substituting for f (x) and for the limits, we have 


apy = th [sind,, (x, + w,)- sind, 7) 
Similarly, 

b,.. = : [cosi\,, x, - Vn OC ea on Po I 

KVP ee rtf, REREMEN Nessa cds 


The overall current density sheet consists of m pulses separated from 


each other by t,. Then 
m m 00 
Nr i 
h(x) = m ho(x) = = : eS (ay, CosA.x + by sind x) 
r Wipe wy Vv Vv Vv Vv 
k=1 ee Se ey at 


(2-2) 


Assuming that the stator conductors are distributed symmetrically about 


the centre line of phase A as shown in Fig. 2-4 


m mM 
br. Ny, ins i hh Wee 
pee ek) keep We sa es 


and a, = 
aie eye Wie ee Fe 
as 
Choosing the origin such that x, = - >> then 
2 m 
Fook Sins (\ Wel 2) 
= hats 2 Ny, cos [A,(k-1)1,] = 


Tr QWy Ay We/2)) gad 


fy ee 1.8 (2-3) 
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Ty 
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‘ a 
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(SS) 


Sunda yileststemaye basudiztelb ors etozoubios 107839@ adi jad goloweaA 
SaS .giZ ot mwoie es A odndg 20 SmEE exsueo ed? 
fat ¥ ys 2 7 ‘ z 
S dal = [ax “os al 2? 2” fey ué 


w 
noid ~~ - = (x ded3 fobs okgtyo of) gatecod. 


= [,rGl-d) A) 809 yy z a | 
Pee | 
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Similarly, for phase B 


m 
iS ine (A. we/Z) 
2 1b Ss is 
Cui ceie ree ons 12) 2 TSG ion te 


- ne (2-4) 
¥37 Yo 


and Kae Sb due to the symmetry of the winding. 


Rotor currents are restricted to the rotor sleeve, the rotor core 


being laminated. In general, the rotor current densities are 


ee = > nes cosh,, xt (amperes/meter2) 
vei 
1 Si r * 3g 2 
Jon by SinA = (amperes/meter“) (2-5) 
vei 


and 


16 


(2-6) 


foe} 

; ed r r ees 

Je = 2 a, cosa, x LaDy, sinh, x 
y= 


FIG. 2-6 STATOR WINDING (one phase) TYPE III 
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| foo aodor ad? (svbale 40307 adi b3 bagontaesx ora eindz7u9 toxof 
ore aglttensb ary xo301 sf2 Taxaadg at shessabeind sete 


(*rsdsc\esx5qms) yaeno Zz. 7 
Ss as 
Sag 
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boa 
(a—S) 42 fate 70+. 7x jke0o 7 a ran es a 
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2.3. Field Equations 
Electromagnetic fields in the stator winding, air gap and the 


rotor sleeve are derived from Maxwell's equations 


vo pe= 0 Op) 

V->D=o (2-8) 

Teen eSs eee (2-9) 
ot 

Vx B= - a (2-10) 

f= soon (2-11) 

B=ud (2-12) 


It is advantageous to use both the vector potential A and scalar 


potential ¢. 


We have 


ot (2213) 


Since the divergence of the curl of any vector vanishes, we can 


satisfy the above by setting 


B =e A (2-14) 
Then = A 
= OB cea re) "k: TT — OA T = a 


“ Byyttenee 


(Of-S) & - =H x ¥ rs 
(11S) fost 7 : 
(S£-8) hu+@ 


teisse bas A lelinsteq soisay sd2 diod sev oF anosuninevbs si ar 
& isttmeteg 
vai aN 


o.. owe 


| 
O=8+V bas 

mpo 8w ,eeiistney tojosv yos 4 Layo 53 io sonegiavih ona si a 7 
gnitise yd eveds sia ytekzee 


(ott) | Eeve® io. 


5 
an 
a 
7 


(eins) om ta +: ot) ee 


a" ths 
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Since the curl of the gradient of any function vanishes, the above can 


be satisfied by setting 


= PN 

E+ yee V oO (2-16) 
therefore 

atatetoAur 

E = aL V o (2-17) 


Substituting (2-14) and (2-17) into (2-9) gives 


9 3A See 


€ ee ovo Veo) P= ty) x Bea 
3 oA = = 
or ue nye A Tae V¢) =pVxH-udJ (2-18) 


To simplify the above, the Lorentz condition 


v7A- pe —sye- ud (2-19) 


For quasi-static conditions, displacement current is negligible and 


equations (2-19) become 


vAze-udJ (2-20) 
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In the quasi-static system the energy that can be attributed to the 
electromagnetic fields is stored in the magnetic field (18). Magnetic 


stored energy in the machine is 


Wn = 5 f B+ H dv (2-21) 
volume 


The stator currents are z-directed so that rotor currents are also z- 


directed. Equations 2-20 become 


v2 Ne =-u ae in the stator conductor (2-22) 
- wu pS in the rotor conductor (2-23) 
and 0 elsewhere (2-24) 


If A, is expressed as a series 
INS SS ALY (2-25) 
vel 


Each component of the series can be found independently due to the 


linearity of the system, therefore 
2 = y2 4 2 % 
ie hemi es ay, (2255) 


The tubular motor has three regions to be considered (Fig. 2-2) 


IV 


(1) the stator conductor Y3 ey Y5 
(2) the air gap ye ayes yy 
(3) the rotor conductor va = y =.0 


(Est) --t6asubngd 197829 St ny * Ea ~ = oh 


(£5-$) :o3¢ubn0o togou oft mk Uy sn 


(RES) Tertupels 0 ; tre ~ _ 


natrae 9 di buednuee Ob gh et 
Sd 


(¢$-S) vg = =A 
ey 
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| a arn. ae 


(St yglt) bersbiacdy of 02 pnotavt vate Ged sosom wilouiee ont 

= pu (+ voilawbmds) abseil 42) 
ges cA gy ang Thue ES) 
02% = Wh yodsubnos ator sit (b) 


th : | 6 Ae 
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Laplace's equation is solved in Appendix A in detail. The following 


expressions are the field equations of the machine. 


(i) Stator conductor 


rm 
tl 


us A, [(K, sinhd,y + K, coshh,y) cosi\, + 


+ (Kk, sinhdyy + K, coshi\,y) sind, x] (2-27) 


jes] 
i] 


y WA [(K, coshiyy 7 Ke, sinhd,y) sinhi,)x - 


- (K, coshi,y + K, sinhdy) cos\,x + 


+ 


Ge sink,x - bo cosi,,x) ] (2-28) 


hy 


(ii) Air gap 


a: Ces sinh y + Ks cosh), y) cos\,,x + 
' 
+ (K, sinhayy + K,, coshi, y) sind, x] (2-29) 
' 2 . 
By = Hay [(K, coshi,y + K, sinhd,y) sind,x - 


™ (K, coshhy + i sinh},y) cos\,x] (2-30) 


(iii) Rotor conductor 


WwW Z Ww 
poe AK, sinhi,y + K, coshi, y) cosi, x + 


Ww Ww , 
+ (K, sinhiy + ie coshi, y) sind,,x] (2-31) 


+ xpheod (xyfitnte 2 +yAdlaar Al) = 


ang tbh (22) 
+ x, ke03 (yAfgoa ts + yArioke pois a, i 


(@f-£) [x Anke (y4h809 2 + y Arik @ + 


~ xybita (yyhiinle (2 + qhrleoo (xy pyhu oe : 


(G€~<) [x Aeoo (vy Atioke yt + ¥,Ateon 3) - 


tmhuor (x Adaoo 2 + ygkeinka TO] 4 = 2B 
| | ‘ r 78 
Ges Exsnte Ce, (t+ ggktiate Se 


vs 
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uw 
BL = HAY CK, coshi\ y + K, sinh\ y) sini)x - 


- (K, coshh,y + Ke sinhi,,y) cosh. x + 


1 ; 
nee (C, sind,x - C, cosh, x) ] (2-32) 
) 
in which 
- ar r a r 
C, = (a,, sinh, x, + by, cosh, X, ) eas 
co (ane cosh, xX," = he sind, ) 


In equations (2-27) to (2-33) subscript is omitted for convenience 


Vv 


for the constants K (i.e. Ky is Ky: 


2.4. Boundary Conditions 


For the given model we have to consider two different interfaces: 
conductor air and conductor iron. We can use the same general approach 
(Fig. 2-7), the result of which has to be modified for the conductor 


iron interface. 


IRON 


CONDUCTOR 


FIG. 2-7 BOUNDARY CONDITIONS 
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Applying Maxwell's equations, we have 
fi: ds=f M+3)- aa 
which becomes 
Ay Ay Ay Ay 


TCM crore Ee ee tet ies iat ee Eves 
xX) ¥3 2 Mig2 Xo : ae 2 Ps 2 ae taste 


If the current density is considered to be finite, then as Ay + 0 


we have 
Bow Axa Ax = 0 
7 


or H =H 
x, x, 


» continuous 


for the conductor-—air interface. 


For the conductor-iron interface, H 0, therefore Hy =H = 
2 


x iron 


It can be shown (16) that the other boundary conditions are 


B = B for the conductor-air interface and B = B for the 
v1 2 bby. ¥3 


conductor-iron interface. 


Applying the boundary conditions to the field equations, the K 


coefficients are evaluated as follows: 
(i) For y-="0, <H. =*0 


te ALK, sinhiy + K, cosh Ay) cosAx + 


+ (K, sinh)y + Ky coshiy) sindx] = 0 


ee ee eee tego ae) 
0 + yh a5 ost otal? aia emma) 


oe gt ral ; 
evounbrtos tll 7 a 


r ' 
sontssant sia soso, ae 
Oe yh = git stotered (0 = 5.48  sastzeint i a aed 


s78 enotitbaos yxabauod rerio ads Jed2 (81) mode 8d mao 37 


oe@) 
aft rot ev - «é bos Sostrevnt its-soaoubnos sit 102 é ; 
.sonizotat a a 


‘avollo? ae baisuleve ox esnataktisog 


© ole gt 0% ear ap 
+ #hans Gh faoo, Ti adite She ge 
O actniite. (eles ela ‘ae 


(ii) 


dole bs) 


gives 
K3 cosAx + Ky, sindx = 0 


which is satisfied for 


H, and B,, continuous, y = Y) 


y. 
Equating (2-31) and (2-29) 


(K, ~ K}) sinhhy, = K3 coshhy, 


1 . 1 
- K,) sinhhy, = K coshhy, 


4 


B,, is continuous, therefore 


' . 
(K, - K,) coshdy, - K, sinhiy, 


- K,) cosh Shee Ky, sinh y, 


H, and BY continuous 


' 
(K, - K,) sinhiy, + (K, ~ K) 
(K, - K,) sinhay, + (K, - K,) 
(K! - K,) coshay, + (K3 - Kg) 


! 
ae) COSI Ay t (K,, 2 K,) 


C 
men 
ye 


r 


coshiy, : 


coshhy, 


sinhhy, 


sinhhy, 
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(2-34) 


(2-35) 


(2-36) 


(2-37) 
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, 13s H, = 0 


I 
oO 


Ky sinhhy, + K. coshiy, ~ 


K, sinhiy, + K, coshiy, = 0 (2-38) 


From equations (2-34) to (2-38) we have 


C 


gest a a 
K, = 5 sinhdy, cothiy, = ae sinhay, cothay. 
ee i See re: hh 
2 = z2 Sinhiy, cothiy, - 12 sinhay, cothay, 
Poets gies Sear eh 
a x2 pera 
igh cas 
AE a ena 
ae ae 
K, = Ky + a coshhy, 
Se ce ha 
Ko = KR, 2 coshiy, 
oe a 
3 x2 1 
ees inh) 
ceeaeaee oa. pa 
Bt eens ee coen 
KR, = K, ae coshhy,, +2 coshay, 
"= K rips h? = £2 cosh 
Be neg ie ee? yy 
eo 
Ky = 
Kk, = 0 (2-39) 
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The set of equations (2-39) lists the K coefficients for the 


field equations (2-27) to (2-33). 
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The parameters of the machine must be found in order to determine the 


equations of motion. First, equivalent inductances of the machine are 


evaluated from the energies stored in electric and magnetic fields. 


zal Magnetic Stored Energy 


Energy stored in electric fields (4) is 


Under quasi-static conditions this is negligible. 


Energy stored in magnetic fields is given by (2-21) 


Mapeciecs = | abet Hedy 
Vv Ply 
eg ta eee 
él 2 il g a 
or Wy = oi J Pa ee = oa J J J (By” + By”) dv 
and Wn = W, + W, + W, 


where Wy > W, and W. are the energies stored within the stator 


winding, air gap and the rotor winding respectively. 


The details of the calculations are given in Appendix B. 


Equation (3-2) becomes 


(3-1) 


(3=2) 
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i 
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ay 


+ 


ea Ee = 
{sinh Yo coth\y, ony ve vo 


= sinh y, cosha 


2 
v¥5] + x7 Hay” Beet Die ea). COS wana 


s reas tc ; “i ¢ 
+ Qn a, a.) by ) sink, xX, ] 


[sinhAy, coshiyy, - sinhdyy, sinhd,y, cothAyy,]} 


(3=3) 


Representing the hyperbolic functions by the first two terms 
of the corresponding series and dropping the fifth and higher order 


terms of 4,,, we have 


- 2 g 2 You 
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Sar Se r Ss fe 
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The magnetic stored energy in the above equation is a function 
of several variables among which are rotor conductor, air gap and 


stator conductor thickness. 


3.2. Inductances and Resistances 

The actual machine is replaced by a coupled system having two 
windings at right angles (1,5) on the phase axes. The rotor and 
stator windings carry currents i,", ivf iS and i,§ respectively, 
the subscripts referring to the winding (phase a, phase b) and the 
superscripts referring to the coil set (stator - s, and rotor - r). 


The magnetic stored energy is then 
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rr r r re P 18 r rr : r 2 
? Laby (i, thy bons Lyay (in, 15.) Daae Loby (iby Nie 
sr s r rs cet Ss sr Ss r 
Penn tayet ay) it) bare (ies tee rt Lele dat 
rs r Ss sr ° Ss r rs : r Ss 
Bs Li ay (i, aay he Poan (ip) tay ae Lapy Cee 1b.) es 
co Leby (ip, +by eee Lyby (iby iby )] 

(3-5) 
where Leu and ee are the self-inductances and inate and ter are the 
mutual inductances. 

Relating currents in (3-5) to the actual current densities, 
the inductances are determined by equating (3-5) and (3-4). 
The rotor coil is assumed to have the same distribution and 
number of turns as the stator winding (1). 
S s Bere 3 r = 
Bj) aac oe) Be Y, a" cos 4 x (3-6) 
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Rewriting (2-3) and (2-4) for the stator current densities 


Ss 
s) =. ay s 
a) ~ = 1, 
¥3 Yo 
x 8 
by, ae iy (3-8) 


Using (3-7) and (3-8), the inductances of the equivalent machine 


are determined by equating (3-4) and (3-5). Then 
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The resistances of both the stator windings are assumed equal, 


k 
P 
R= R,§ aR eee Ree ior Ny R (of the conductor) 


j=l 
(3-10) 
and are associated with the first space harmonic (1), therefore 
ae a - 
te Ruy Raa for, Wee a 
= 0 fone yer ok (3-11) 


The equivalent resistances of the rotor are found by equating 
rotor losses in the actual and equivalent machines. The total 


instantaneous rotor power dissipation in a machine is 


r 2 2 2 
Pa= & f ff | ah + aj0 45,7 | ax® ay ae (3-12) 


which for our model becomes 
» he i 
p 2 2 
Pg = = per f i ide 4 Jt ) dx*® dy (3-13) 
0 O 


Integrating the above equation, the equivalent rotor resistances are 


R Ys Rea sz ee eo ool (3-14) 


The equivalent inductances and resistances derived in this 
chapter are coil variables dependent upon the relative positions of 


the rotor and stator axes. 
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CHAPTER IV 
DEVELOPMENT OF FORCE EQUATIONS 


eral GA Equations of Motion 


Equations of motion are best expressed in matrix notation, which 
is used throughout the chapter. As done in (1), it is convenient to 
subscript the matrix rather than the individual harmonic variables in 
the matrix. 


The volt-amp equations are given by 


where 


(4-2) 
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The electromagnetic force on the 
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eave 


rotor (1) is 
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Assuming a linear magnetic circuit 
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and 
0 0 | -sindx’ -cosdxt 
of 0 0 cosAx® -sin)xt 
|S REE ek, 
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v 
(4-9) 


4.2. dq Transformation 
Equations (4-1) and (4-7) describe the behaviour of the machine 
under all operating conditions, however, they are non-linear. 
Applying a linear dq transformation, the space dependance of 
the equivalent inductances and resistances is removed. The trans- 


formation is power invariant (5) and is defined as 
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The following equation holds for the transformation matrix: 


a = a (4-13) 


18 


Equation (4-10) implies that the set of rotor currents i 


7 6% 
and iny 


in the rotor windings moving with the stator is transformed 
to the currents shee and i which are stationary with respect to the 
stator reference frame. It means that the air gap fields produced 
by currents ice and ie are the same as the fields produced by the 


currents be and i,* flowing in the fictitious windings stationary 


to the stator. 


We have 
v ie (4-14) 
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Q . L a = 
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Equation (4-7) becomes 


ray = = Ca, eel re F aqv 
* + ia (a, ae a,) © aqv 
ie = Le AN a ans = dqv (4-20) 
or 
aS alee gcd S Pog ae. (4-21) 


Equations (4-19) and (4-21) constitute the set of equations of motion 


of the machine. 


4.3. Steady State Relationships 


The fundamental slip is defined as 
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For steady state 


p = ju° 


and the volt-amp equations are 
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Considering the steady state and balanced conditions, the 
Ss e Ss e 
stator currents i. | and ip) are in quadrature. 
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and 


Ce rte ae) (4-26) 


Similarly for the currents we have, dropping the subscripts, 
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The force equations become 
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where 


Ss = =r VS, forward harmonic slip 


n 
i] 
| 
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S 

n 


backward harmonic slip (4-32) 


4.4. Equivalent Circuit 
Equations (4-31) can be represented by the equivalent circuit 


shown in Fig. 4-1. 
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FIG. 4-1 EQUIVALENT CIRCUIT - STEADY STATE 
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CHAPTE RI SV 


DIGITAL STUDIES OF THE MOTOR 


The entry and exit edge effect of the machine are neglected, there- 
fore calculations of the machine performance are based on the straight- 


forward use of equations (3-9), (3-11), (3-14), (4-29) and (4-31). 


5.1. Force Equations 


Harmonic forces can be derived with the help of the equivalent 
circuit or directly from equations (4-31) and (4-29), giving the same 
result. 


The total impedance of the equivalent circuit is 
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The kth harmonic rotor current is given by 
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and the kt) harmonic force by 
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The total force on the rotor is obtained by summing the 


harmonic forces (5-5) 
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In general, some of the spatial harmonics of the current 
densities are negative, creating backward travelling magnetic fields, 


therefore the corresponding harmonic forces in (5-6) are reversed. 
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ee Efficiency 


The performance of the machine was computed on a digital com- 
puter, using the APL program language. The machine data are listed 


in Appendix C-l. 


Efficiencies are calculated, using the following equations: 
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power output + losses (5-7) 


Hysteresis and iron losses as well as mechanical losses of the 
motor are not included. 
Winding distribution factors for different types of stator 


windings, digital programs and inductances are listed in Appencix C. 
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CHAPTER VI 


RESULTS AND DISCUSSIONS 


6.1. Outline 

The steady state output characteristics of a tubular motor were 
computed, using the program in Table 2, Appendix C. The program holds 
for an infinitely long machine. As a tubular motor would operate at 
higher slips, the entry and exit effects are negligible (20) as both 
the forward and backward travelling waves caused by the reflection at 
the edges decay rapidly. 

The effect of different winding distribution factors, air gap 
thickness and rotor sleeve thickness on machine characteristics is 
studied in this chapter. 

The calculations are done for three stators, each having 
different winding distribution. They are type I, type II and type III 
stator windings, as described in Appendix C-2. Thickness of the 
winding is given by the value of ly3 - y>| in centimetres. A set of 
different rotors is considered, each having different thickness of 
the copper rotor sleeve. Another similar set of rotors has aluminum 
sleeves. The thickness of the rotor sleeve is described by the value 
of Y in centimetres. The clearance, i.e. the air gap between the 
rotor sleeve and the stator winding is ly, aay and is also in centi- 
metres. 

The machine data are listed in Appendix C. The temperature of 


both the stator and rotor conductors is assumed to be 75°C. 
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6.2. Effect of Winding Distribution Factors 

In Figs. 6-1, 6-2 and 6-3 steady state output forces, effi- 
ciencies and stator currents are plotted. 

Fig. 6-1 compares output forces for three different stators and 
the rotor with the sleeve 0.2 cm thick. The characteristics for the 
copper rotor sleeve show an increase in output forces as the stator 
windings type II and type III are considered. The type III curve 
reaches a maximum for slip s > 1. The aluminum rotor sleeve has 
higher rotor resistances, therefore the force maxima occur at slips 
Sel; 

Fig. 6-2 shows the effect of stator distribution factors on 
efficiencies. For the same type of stator winding, the efficiency is 
higher for the machine with the copper rotor sleeve over one with the 
aluminum sleeve as the rotor losses are smaller for the copper sleeve. 
The efficiency of the type III machine is higher than the efficiency 
of the type II machine, which is higher than the type I machine effi- 
ciency, as the distribution of magnetic fields in the air gap is im 
proved. 

Fig. 6-3 shows the stator currents of the machine. The highest 
currents are for the type I machine, the lowest for the type III 
machine. 

The performance of the machine was considerably improved by the 
type II stator winding and was the best for the type III stator 
winding, as the direct result of the improved distribution of the 
magnetic fields within the air gap, the distribution being more 


sinusoidal. 
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The following data were used: 


ee 0.2 "cm, yo 0.3 cm 


Figs. 6-4, 6-5 and 6-6 are similar to those considered in the above 
example. The thickness of the rotor sleeve was increased, thus the 
air gap y3 was also increased. The rotor resistances decreased, 
therefore the force maxima occur at lower slip. The machine per- 


formance is as described above. The data used are 


‘btn aks Ozosen. Yo = 0.4 cm 


Fig. 6-7 compares the first harmonic force with the total force, the 
difference is less than 3% for slips higher than 0.2. This is for 
the machine with the type I stator winding distribution. The over- 
all effect of the higher harmonics is small and the total force is 
smaller than the first harmonic force. The effect of higher harmonic 
forces for the type II and type III windings was much smaller than 


for the type I winding. 
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6.3. Effect of Air Gap Thickness 


In the next example only the winding type I is taken into con- 
sideration, as it is the simplest to make. The air gap thickness 


ibe which is actually the clearance between the stator and 


ie 
rotor conductors, was assigned different values and the program in 
Table 2, Appendix C, was used for computations. 

Fig. 6-8 shows the effect of variations of the air gap on the 
output forces both for copper and aluminum rotor sleeves. An increase 
of the air gap width causes a decrease of the output forces. This 
is due to the increase of the leakage of the magnetic fields within 
the air gap. The curves corresponding to the aluminum rotor sleeve 
have the peak forces near standstill due to higher resistances. 

The efficiency characteristics for the tubular motor with the 
copper rotor sleeve are shown in Fig. 6-9 and with the aluminum sleeve 
in Fig. 6-10. The decrease of the efficiency is considerable except 
for speeds close to standstill. 

The stator currents are plotted in Figs. 6-11 and 6-12. The in- 
crease in the air gap width decreases the inductances as can be 
seen from equations (3-13). This means that the total impedance of the 
machine is smaller, therefore the stator currents are larger as the 
terminal voltage is constant. Furthermore, the leakage inductances are 
increased as the air gap width is increased. The magnetizing induc- 
tances are decreased, therefore requiring higher magnetizing currents. 

The machine performance is considerably affected by the increase 


of the air gap, as both the output forces and efficiencies are de- 


creased and the stator currents are increased. 
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6.4, Effect of Rotor Thickness 

Next the set of tubular motors, having the constant air gap Y3 
and variable rotor thickness Yi> is considered. Fig. 6-13 and Fig. 
6-14 show the output forces for copper and aluminum sleeves respec- 
tively. For small slips the forces are increased, for higher slips 
they are smaller. The efficiency characteristics (Fig. 6-15 and Fig. 
6-16) and stator currents (Fig. 6-17 and Fig. 6-18) are affected in a 
similar way as the forces. 

The increase of forces and efficiencies is the result of in- 
creased rotor conductivity. As the slip increases, the effect of 
increased leakage inductances is larger and the forces and efficiencies 
decrease. This is more apparent with the copper rotor sleeve than 
with the aluminum rotor sleeve with higher resistivity. 

The tubular motor, employing rotors with different thickness 
of the sleeves, has better force and efficiency characteristics at 
higher speeds of the rotor. As the direct result of increased rotor 
leakage inductances for small slips, the increase of the rotor sleeve 
thickness causes the decrease of the output forces and efficiencies. 
For small slips the machine performs better, as the effect of higher 
conductivities is bigger than the effect of increased rotor leakage 


inductances. 


» Weare ee ” 
en . : Sal 


7 oe ’ , + - . 
_ 7 : a / ’ » 7 n 
= a7 - — 2 - 
. ' ‘apsmla bet rodot OITA Mow 
ich a 


i 7 : - : : | 7 

eX Ga ths Jonteuos sit ankved ,axotom x6 Ludu3 to Jee os 3 
; a 2 ee re = eek 
gr? bem €2-0 ght .bazabtanos st . /¥ azacdoinys t0d0% o Like bea bas 


= 


: < 
-seqae? seveele muntnate bos. 3» 9qnoo “ot noniot tuqzuo ond wore 


agiba vertgid zot ,bsdaszoot wis ssoTot il) eqiie Iflawa vod “levis , 


efi baa CL-O .giT) estiaivosoureio yonsiotitoe eM celia 63a yoda — 
7 


n oi bagosils ste (81-4 .9i bone Vi-0 Jalt) etoorrn Todese bem (05-3 


-2eo7t0t od? ew you talinis 
~a) %o sfgeer Sa oY astonslritis bes <s21Tc? Jo, sansronk eft 


ts gootte sid ,eseesyonr qita and #A »vitvidovbaoa 10305 beaseTs 


+ 
i 


esisteioltis bas essto? of? bis wogyel «! asonesodbn! sgetest beasstoo 
aatiy oveste vote? seqqo> std Wotv moveqqs stow at aps .sanotosb 


iivisateos ssngti d4fw ovesle oder muninults od3 datw 


guendodds sastel2tb dytw aroyor gatvolqans ,r0J0m usludus giiT 


ih pokiulreieareds yous larYio bas sore? zedied de 39 peala atta De 


fore) bsaasdont Yo tikes Jotatb oid eA =. t0tDy gf Jo Rbaeqe Ts itd 


. 4 ~ 
veel rtodor oy to sandtort sis agile I lemte rol sssdassuont mods b 
- ¥ 
seetoes i stits bab asorct juqivo sty to oa bb sty pazval €R oy 


apiigiyi © Yastie of3 as .35)394 emroiupd eotdorm add agile (igera- 
ne 


gupteel wciot BaeseYoni to Joddin alle oua ¥ agqtd fi on 


a a 


56 


dil'IS 


p 


ut ¢ 


WU G*y fy FG°E FE FG°z 


SAZ2TS 


nO 


SSHNMOIHL YOLOY AO LOgdda - 


OF 8°0O oEG 


Sd0u0d IAdLNO €T-9 “DIA 


7°0 G0 


OS 


OOT 


OST 


002 


[SNOLMAN] dOuod 


Saou0d INdINO ¥YI-9 “Old 


wu ¢ = A 


[SNOLMAN] dOuOd 


um ¢*y th cre te SG*c = A 


eAs2eTS TV 


58 


pepnToUT Jou sesso, TeoTueyoou 


mu ¢ OK 


I 


wml ¢*h SC°E 2O°Z A 


dAsEeTS NJ 


AONHIOIdde GT-9 “OIA 


[%] AONAIOLAAg 


Bo 


pepnpTout Jou sessoT TeoTueyoow 
wu G¢ = we 

6 ° € ry I 
Mage eae Gr. ame G 


SAZ2XTS TV 


AONSIOIdda 91-9 “OTA 


[%] AONAIOLAIa 


avesie IA 


“gi 2.8 12,8 32.8 


: rom = 
_ 


c- _ bebuiont Som doeeol Taotasdoom _ 


Ss 
7 . Za 
a c.f —— 
oo es ee ee 
a 7 — ™ ~~“ 
- > -—~ ee — ~ \ 
> a ~ 
_ — > <a Sas SA 
eS i ~ 
5 Ee esi 
qr12 0.1 &.0 3,0 6.0 $.0 


YOMAIOTTIA Of-0 .OT7 


60 


wu 


UW Gy 56°C +O °7 


BAd2TS 


nO 


OF 


8°0 


SLNAWUND YOLVLS LT-9 


9°0 


‘Old 


OL 


ct 


[Sdtv] LNdwano 


SLNAMMND YOLVLS 8-9 “Old 


dI 1S O°T 8°0 9°0 7.0 c 0 
aa $$ $f $$$ —___j— 
ag 
Sf ee Se 
a ot 
Wu G = eK 


PAVOTS TV 


[SdWv¥] LNaaano 


62 


6.5. Constant Clearance 

Next, for practical purposes, the air gap clearance ly, raph ie 
between the windings is kept constant. The effect of the rotor con- 
ductor thickness is considered for the type I stator windings. This 
means that an increase of the sleeve thickness causes a corresponding 
increase of the total air gap thickness, thus the effect of both on 
the machine performance is combined. The air gap clearance is 0.1 cm. 

The output forces are plotted in Fig. 6-19 for copper sleeves 
and in Fig. 6-20 for aluminum sleeves. An increase of the sleeve 
thickness causes an increase of output forces for small slips and a 
decrease of the forces for higher slips. The decrease is caused by 
the increase of the leakage inductances. 

The efficiency curves in Fig. 6-21 and Fig. 6-22 show some in- 
crease in efficiency for small slips and decrease for higher slips, 
being the result of a combination of rotor sleeve resistances and 
leakage inductance. Fig. 6-23 and Fig. 6-24 show an increase of stator 
currents as the sleeve thickness is increased and consequently the total 
impedance of the machine is decreased, as the increase of leakage in- 
ductances has a smaller effect than the increase of the conductivity. 

The machine performance is affected for the higher slips, being 
decreased especially for the machine with the copper rotor sleeve. 
Therefore, an increase of the sleeve thickness, while the clearance 


is kept constant, results in a machine with inferior output cha- 


racteristics. 
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6.6. Conclusions 

A tubular motor is a linear induction machine which is very 
simple to make. This holds especially for a machine with the type I 
stator winding distribution. Motors with type II and type III 
winding distributions are more difficult to make, but the resulting 
increase of efficiencies is considerable. This increase is approx- 
imately 40% for type II and 100% for type III stator windings over 
type I windings in our example. It is also interesting to notice the 
drop in the stator currents for type II and type III windings, caused 
by the increased inductances and rotor resistances, i.e. the total 
impedance of the machine is increased. 

When designing a tubular motor, one should take the cost of 
winding and the increase of efficiencies into consideration. 

Next the machine with a constant rotor sleeve thickness was 
considered. The air gap thickness variations showed that an increase 
of the air gap causes an increase of leakage inductances, thus 
deteriorating the machine performance. In other words, a designer of 
a tubular motor should keep the air gap as small as possible. The 
motor with fewer layers of wires in the stator winding and with the 


larger repeatable section 1, would have higher efficiencies. 


Yr 
When the air gap thickness is kept constant and the rotor sleeve 

thickness is assigned different values, while the clearance is changed 

accordingly, the machine performance is dependent on the slip. For 


small slips, the increase of rotor conductivity has a bigger effect 


than the increase of the leakage inductances and the machine performs 
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better. For higher slips, however, the increased leakage inductances 
cause decrease of the machine performance. This means that the 
machine having smaller amount of rotor conductive material could per- 
form better than a machine having a thicker rotor sleeve when the air 
gap thickness is concerned. 

When a tubular motor has a specified clearance and a stator 
winding, one can change the output characteristics, using rotor with 
different thicknesses of the rotor sleeves. Now the air gap thickness 
and the rotor sleeve thickness are combined and the machine performs 
better for low slips, as the rotor conductor thickness is increased. 
For higher slips, the effect of increased leakage inductances is 
bigger than the effect of increased rotor conductivity and the machine 
performance is decreased. Therefore, when designing a tubular motor 
for low speeds with the clearance specified, one should design the 
sleeve as thin as possible. There is a limit to this, as comparison 
of Fig. 6-21 and Fig. 6-22 reveals that the increase of rotor resis- 
tances causes the expected decrease of the efficiencies. 

The effect of higher harmonics was found negligible for the 
applied approximation of the hyperbolic functions for the assumed 
stator winding distributions. For different types of windings the 
approximation may not be satisfactory. 

The theory used within the paper is general and may be used 
for LIMs. The inductance equations (3-9) hold for machines 


satisfying the assumptions in Chapter II. 
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The study presented above is applicable to the "infinitely long" 
tubular machines. In practice, this means that the more pole pairs 
the machine has, the more exact the equations are. The assumption of 
an infinitely long machine lies in equations (3-6). To solve the 
entry and exit edge effects, the different boundary conditions would 
have to be found; this in turn would lead to more complicated equa- 
tions. It would be useful to develop the above theory by the inclu- 
sion of the edge effects, as the machine performance could be pre- 


dicted more exactly, especially for low slips of the rotor. 
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APPENDIX A 


Solution of Laplace's Equation 
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There are 3 regions to be considered for the tubular motor 


(Fig. 2-2): 
(1) the stator conductor Y3 2y2 - 
(2) the air gap Yo 2 y2 yy 
(3) the rotor conductor yi, =y 70 


(i) The Laplace's equation (2-24) in the air gap 


927A 927A 
VAA = Roo + a > y2 A,y = @) Vv = L a 3 e 


oy =j 


has the solution of the form 


A = A(x,y) = X(x) Y(y) 

Then a*A7 " 
ee (x) 
927A 


= X(x) Y (y) 


and (2-24) becomes 
YX) xy! o= 0 


or " 


x" Y 
Pn) a 


(A-2) 


(A-3) 


(A-4) 


O5¢5 -% 


qag tie sid ut (S-S) amtteups s'soatqnl od? 


2 = AS <. s ‘ = Sy 
ere) Aa Pe uv o ysh Y = arg + PS xh 
(I-A) 7 | ~ 
sro? ont 20 notieloe od? end 
{2-A) «AQT Get = Coda = A 
gs cont 
| al * (oY “ier 
‘ A*6 
(¥) ¥ (x= See = 
(Ena) 


0 = "YS + ""Y 
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Since x and y are independent variables, they do not depend on 


each other and (2-24) is true only if each term in (A-4) is the same 


constant, therefore 


X ni=iac X 

W 

Veh =) = iChy, 
Let c=- 2, then 

Ww 

ee he 

xX 

W 
and ae he 

Y 
Therefore 

Ww 

X. = 2. X =.0 

Ww Q) 
and Y (2 =. A> -Y 


Solving the equations, we get 


A = X(x) Y(y) = (By sinhay + B, coshAy) (C sindx + 


+ D cosdAx) 


Then 


Zz 


+ Ky, sinhidy) sindx] 


= a a 2s JA, . jA 
Since Be eS oe ae ee ee ean xt 
oA 
we Boe ee aa 


i 
AJ Fu (kK, coshiy + K, sinhay) cosAx + (K,, coshay + 


(A=5) 


+ xfhnie 3D) (yAdeos af + ¢khiate A) = (y)¥ @)X =A 
kaos G + 


| + yineod 2» + -xheg> (ehkdake a + wkdeoo > y= Po 
_ A) leknke (ysilnte (H+ 


(ii) 


ial 


and u HH. =-—=B (A-6) 


Therefore, from (A-5) 


mo 
tl 


x = AL(K, sinhdy + K, coshy) coshx + (K, sinhay + 


+ K, coshiy) sindx] oe) 


ee] 
i 


y WALK) coshhy + K, sinhdy) sindx - (K, coshay + 


+ K, sinhay) cosAx] (A-8) 


Region 1 - the stator conductor within which (2-22) holds: 


ghee “OA, 
Sout aD 
oy 


: = - u(a® cosikx® + b® sindxS) 
x 


and x® = x, y® = y, z® = z, so we have 

397A 927A 

Se + = = - u(a® cosix + b® sinAx) (A-9) 
ox Oy 


which is satisfied by 


A, = ul(K, coshiy 4 K, sinhidy) cosirx + (K, coshiy + 


zr sinhdy) sindx hy (a® cosax + b& sindx)] (A-10) 
Therefore 


H, = A[(K, sinhay + K, coshiy) cosdx + (K, sinhdy + 


+o, coshiy) sin\x] (A-11) 
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By = WAT(K, coshhy + K, sinhay) sindx - (K, coshAy + 


+ K, sinhdy) cos\x ie (a® sindx - b® cosix)] (A-12) 


Region 3 - the rotor conductor in which (2-23) holds. 


Referring the rotor currents to the stator reference frame 


a 
bmi J 
zt=z 
x x Xo x > X_ 


where x,* is the rotor origin referred to the stator and 


Kee = x™t for rotor motion gives the Galilean transformation 
(13.5°73)-. 

We have 

eam ae 

— = + —+ = - u(a™ cosix™ + b? sinix‘) = 

ox oy 


= lat cos) (x,-)x,7) + b* sink(x -7x,*) ) = 


- ulcosrx(a* cosAx,” Eth replay 4809) Pa 
+ sindx(at sindx,* eb ~ cosAx,*)] 
The solution is of the form 
' W 
Neue ul(K, coshiy + Ke sinhdy) cosax + (K, coshiy + 


5 Ke sinhidy) sindx] fa [(a cosax,* - bt sin\x,*) cosax + 


+ (a® sinix,* Sus cosAx,*) sinAx] (A-13) 
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Therefore 


Ww 
By = - Au[(K, coshdy + Kj, sinhdy) cosdx - 


(Ky coshiay + ro sinhiy) sindx + 


1 
bub ¢at sin\x,* + pt COSAX,") COSAX - 


re cosAx,* abe sinhx,*) sindx) ] (A-14) 


Hy, = LR, sinhdy + K, coshAy) cosAx + 


= (K, sinhhy + K, coshiy) sindx] (A-15) 


let 


(at sin\x,* dab “Cos\ xX a) 


i} 
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i] 
(2) 


(at coshx,~ - pt sinhx,”) (A-16) 


Substituting (A-16) into (A-14) 
i] 
By = vrL(K; coshhy + Ke sinhidy) sindx - 


- (K, coshiy + K, sinhiy) cosdx + 


++ (C. sindx - C, cosix)] (A-17) 
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APPENDIX B 


Stored Energy Calculations 


(i) Evaluation of (3-2) gives 


v= aa » oe he taker eae Gina y ee ecetilay 
vel 


- sinhAy, coshd,y,) + 2(K,\K ) (sinh?A OSs - sinh?) y,) + 


1Vv 3y 7 Ko yk 4yvVv vo 


2 s : ; 2 s 
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ak 32 52 _ = 
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Lge. Sum G2) vain t= ee eT (B-3) 
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Because of the complexity of the above expressions, first the 
hyperbolic functions were represented by the first terms of the 


corresponding series. If the ratio of air gap thickness y, to 
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wavelength ty is small (17), then coshdy, * 1 and sinh y, = 0.1 and 


sinhay3 ~ dy3, sinhdy, * Ayo, sinhdy, * dy). 


Also sinhhy, s sinhhy, aS sinhdy , 


1A 


1 = coshay, = coshiy, 4 coshiy, * 1 


Thus the set of equations (2-41) becomes 
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Substituting (B-4) into (B-1) to (B-3), we have 
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RuP 1 2 aly 
Y= ae » - Bee et abS — (y, -y,)° + 
v= ws 3 
2 eit 2 
eer re Cea a oe ti ad a ek ta: 
+ 2) (Co ac + s zal 
ne Cees sy) o y, ae a | (B=7) 


Inductances of the equivalent machine were determined from the 
above equations as 
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resulting in Ley, > L > L and negative leakage inductances of 
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both the rotor and stator. Therefore the above approach was not used. 


(ii) Instead, equations (2-40) were substituted into (B-1) to (B-3), 


resulting in the following: 
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foe) 
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v= 
+ sinh*Ay,) (sinhdy, coshiy, - sinhdy, coshiy,) - 
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APPENDIX C 


MACHINE DATA AND CALCULATIONS 


C-1 Machine Data 

The machine can employ different rotors, as they are easy to 
replace. This is a big advantage of tubular motors, as different 
characteristics of the machine can be obtained simply by inserting 
rotors with different thickness of the conducting sleeve and/or 
different sleeve material. 

The idealized machine used for digital analysis has the 


following data: 


Overall diameter of machine 7 cm 

Thickness of stator winding 0.6005 cm 

Stator length 24.2. cm 

Power supply 110V, two phase, 60 Hz 

Stator winding 432 turns per T, per phase 

Wire #18, 7.76 2 per 1000' at 75°C 

Rotor sleeve copper, o = 4.93 x 10’ mhos/metre at 75°C 


aluminum, o = 2.97 x 10’ mhos/metre at 75°C 


Four pole stator TegGedeed CtyMe Fit pe 7.26 m/sec. 
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A simple stator winding such as one having the current density 


distribution shown in Fig. 2-3 will be called type I stator winding. 


The distribution factors are derived from (2-3), giving 


gl 
Ke = = eee cos = 
aN Tr aL 4 
2 


where N is the number of conductors within the shaded area in 


Fig. 2-3. 


We have 


a s 
Kay Roy 


due to the symmetry of the windings. 


(C-1) 


Type II stator winding has the stator current density distri- 


bution as shown in Fig. 2-5. From equation (2-3) we have 
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Winding is symmetric, therefore 


Type III stator windings have current density distributions 


shown in Fig. 2-6. Substituting into (2-3) 


a a 
sin — 
s _ 4N 12 TV TY 
Ke a see {3 + 2 cos ae cos > 
a2 
- cos it - 2 cos SL) 
3 6 (c-3) 


where N is the number of conductors within the shaded area in 


Fig. 2-6. 
We have again 
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The winding distribution factors for three types of stator 
winding distributions are listed in Table 1. The amplitude of the 
first harmonic is the largest and the machine will rum on this pre- 


dominant harmonic. 
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C-3 Digital Programs 

The program in Table 2 is written for inductances derived from 
equation (3-3) for magnetic stored energy. The hyperbolic functions 
are preserved and therefore even higher harmonic inductances are 
exact (for the assumptions in Chapter II). 

Table 3 lists the program based on the equation (3-4) for the 
inductances where approximations were used. They introduce an error 
into the higher harmonic inductances. 

Both programs in Table 2 and Table 3 are used for the type I 
winding distribution. The calculations of type II and type III 
winding distribution machines were done with the programs in Table 4 
and Table 5 respectively, in which approximations of hyperbolic 


functions were used. 
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C-4 Inductances 

A comparison between the two programs showed, as expected, an 
increase of the error as the air gap was increased. For large air 
gaps (i.e. y, = 0.25 cm, y, = 0.5 cm), the error for the first harmonic 
inductances was less than 0.2%, while for the third harmonic it was 
almost 10%. The overall effect on the machine performance was negli- 
gible, the only considerable error being for zero slip as the effect 
of higher harmonics was stressed. 

For machines having very large air gaps or having stator winding 
distribution with considerable higher harmonics or for machines running 
on the higher harmonic, the program with hyperbolic functions must be 
used. 

As the values of the air gap and rotor sleeve were changed, the 
machine parameters changed too, therefore, next is an example of para- 
meters for ea 0.2 cm and Vea 0.3 cm only. The machine has a type I 


winding. 
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TABLE 6 ROTOR RESISTANCES 
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